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Spectral sums of the Dirac-Wilson operator and their relation to the Polyakov loop are thoroughly investi- 
gated. The approach by Gattringer is generalized to mode sums which reconstruct the Polyakov loop locally. 
This opens the possibility to study the mode sum approximation to the Polyakov loop correlator. The approach is 
rederived for the ab initio continuum formulation of Yang-Mills theories, and the convergence of the mode sum 
is studied in detail. The mode sums are then explicitly calculated for the Schwinger model and SU (2) gauge 
theory in a homogeneous background field. Using SU(2) lattice gauge theory, the IR dominated mode sums 
are considered and the mode sum approximation to the static quark antiquark potential is obtained numerically. 
We find a good agreement between the mode sum approximation and the static potential at large distances for 
the confinement and the high temperature plasma phase. 
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I. INTRODUCTION 

Color confinement and spontaneous chiral symmetry break- 
ing are the two most relevant features of Yang-Mills theory 
when the structure of matter under normal conditions is ex- 
plored. Both phenomenons are attributed to the low energy 
sector of Yang-Mills theory, an analytic description of which 
is hardly feasible due to strong couplings between the basic 
degrees of freedom, quarks and gluons. Lattice gauge simula- 
tions found that both so different phenomenons are intimately 
related: in the chiral limit, the critical temperature for decon- 
finement Td and the critical temperature for chiral restoration 
T c coincide yj. This finding is highly nontrivial since the 
expectations had been that two totally different mechanisms 
involving different energy scales were at work for each phe- 
nomenon. Indeed, if quarks which transform under the adjoint 
representation were considered instead of the quark fields of 
standard QCD, largely different values for Td and T c were re- 
ported in 131. 

These results have stirred the hope that a single low en- 
ergy effective degree of freedom is responsible for both phe- 
nomenons. Given that in pure SU (N) gauge theory the long 
distance part of the static quark potential depends on the N- 
ality of the quarks and that the crucial difference between fun- 
damental and adjoint quarks is iV-ality again, center vortices 
appear as a natural candidate for such a degree of freedom: 
these vortices are tightly related to confinement, are sensible 
in the continuum limit and offer an intriguing picture of de- 
confinement at high temperatures iH 0, d, |g] . It was found 
quite recently that the vortices extend their reach to a descrip- 
tion of spontaneous chiral symmetry breaking as well 0]. 

In order to reveal a model independent link between con- 
finement and chiral symmetry breaking, Gattringer proposed 
to reconstruct the Polyakov loop expectation value (P) in 
terms of a particular spectral function 6>jv t of the lattice Dirac 
operator QSQ . While the low lying modes of the Dirac opera- 
tor are directly related to chiral symmetry breaking by virtue 
of the Banks-Casher relation, the Polyakov loop expectation 



value serves as the litmus paper for confinement. Spectral rep- 
resentations of (P) have subsequently been the subject of re- 
cent studies. Bruckmann et al. ]9[] investigated the response 
of the eigenvalues of the staggered Dirac operator to a twist of 
the boundary conditions. They found that the infrared part of 
the spectrum is most sensitive to twists. It was subsequently 
pointed out that the spectral sum Sjv t , originally proposed by 
Gattringer, is dominated by the ultraviolet part of the spec- 
trum since the sum contains large powers of the eigenval- 
ues ifsL l9l ITol [Till I . A sensible continuum limit was caught into 
question. In lUOtl alternative spectral sums were put forward 
which serve as order parameters for confinement and which 
receive their main contributions from the infrared part of the 
spectrum. One of these sums is of particular interest since 
it relates the dressed Polyakov loops to the chiral condensate 
via the celebrated Banks-Casher relation lfl2ll . If one twists 
the gauge field by a complex number z with unit modulus 
(or equivalently twists the boundary conditions by 1/z) and 
picks the coefficient of z k in the spectral sums one obtains the 
dressed Polyakov loops with winding number k ifioll . With 
this method Bilgici et al. lfl3ll connect the eigenvalue den- 
sity at zero and therefore the chiral condensate to the dressed 
Polyakov loops. Numerical results for spectral sums of vari- 
ous lattice Dirac operators with quen ched configurations can 
be found in the recent papers if^L 1 1 OL [T2I [Till an d for ensem- 
bles generated with dynamical fermions in 11111 . We finally 
point out that the spectral approach to the Polyakov line has 
been extended by using eigenmodes of the Laplacian opera- 
tor iTTsfl - There, the spectral sum acts as gauge invariant low 
energy filter which reveals the "classical" texture while quark 
confinement is still active. 

In this paper, we generalize the concept of the mode sum 
approximation to reconstruct the Polyakov loop locally. This 
has the great advantage that now the Polyakov loop correla- 
tor and the static quark antiquark potential can be studied in 
the light of a few low lying modes of the Dirac operator. We 
then point out that the mode sum approach is not solely tied 
to lattice quark operators and present an explicit construction 
of this approach in the ab initio continuum formulation of 
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Yang-Mills theory. The convergence of the mode sums for 
all polynomial functions of the continuum Dirac operator is 
demonstrated for the first time. We then argue that a wide 
class of IR dominated mode sums are in fact proportional to 
the Polyakov loop. This conjecture is fostered by an explicit 
calculation of these mode sums for Schwinger model and for 
SU(2) gauge theory with constant background field strength. 
We then consider SU(2) gauge theory above and below the 
deconfinement temperature by means of lattice gauge simula- 
tions. Most important, we find that a few low lying modes of 
the quark operator are sufficient to reconstruct the static quark 
potential at large quark antiquark distances. This is a gauge 
invariant and model independent signal that the color confine- 
ment mechanism has its fingerprints in the low lying quark 
spectrum. 



H. RECALL OF SPECTRAL SUMS FOR LATTICE 
MODELS 

The j-th power of a Dirac operator T>jj with nearest neigh- 
bor interaction on an Euclidean lattice with N t x iV^ -1 sites 
can be expanded in Wilson loops of length up to j, 

(x\V&\x) = J2 a ^ W ^- (1) 

\C<c\<j 

The value of the coefficient ac T multiplying the holonomy 
Wc x of the loop C x with base x depends on the type of 
fermions under consideration. The expansion (fTJ is used to 
relate the Polyakov loop 

P(x)=txV(x), V(x)=Y[Uo(x ,x) (2) 

with the spectrum of the Dirac operator @]. On a lattice 
with N s 3> Nt the sum on the right hand side of (Q~|) con- 
tains contractable loops and loops winding once or several 
times around the torus in time direction. In order to re- 
late the Polyakov loop to spectral sums of quark eigenmodes 
an interface is inserted into the lattice gauge configuration 
Ufiix), x = (x , x) by 

„. , > _ j z ■ Uo(xo, x) for n — and xo = 
U„{x ,x) otherwise , 

(3) 

where z = e 2lrM is a center element. Contractable Wilson 
loops are invariant when one inserts an interface which is re- 
ferred to as twisting the gauge field. Wilson loops winding 
fc-times acquire a factor z k . It follows that for a twisted gauge 
field the coefficient of z in the series (Q~|i becomes a linear com- 
bination of dressed Polyakov loops passing through x, having 
length < j and winding 1 + kN times around the time direc- 
tion. Here k is an integer and z N = 1. For j = N t there is 
only one such loop, namely the straight Polyakov loop at x. 
After taking the traces over spinor and color-indices the sum 



over x° yields JH] 

N n p 

Kj 

k=i p=i 

\(x)=J2\( x °> x )- ( 5 ) 

The first sum in (0]i is over all center elements z\ , . . . , zn and 
the second over all n p eigenvalues of the Dirac operator. The 
value of the constant n depends on the type of lattice Dirac 
operator under consideration. z Q p (x) is the eigenvalue density 
and z X p the p-th eigenvalue of the Dirac operator z T>u = D^u 
with z-twisted gauge field, 

( z Vu) z ^ p = z \/^ p with 

z 4>p(x + N t , x) = - z ip p {x , x). 

In terms of the normalized eigenmodes the color-blind density 
reads 

%(x) = ^2\^ Pti (x)\ 2 , (7) 

t 

where the sum extends over all eigenfunctions of z T>u with 
fixed energy A p . The densities z g p are gauge invariant scalar 
fields. For the trivial center element z = 1 we often write 
q instead of z g. Averaging the local identity over space 
yields 

s x k p 

This simple formula for the averaged loop has been investi- 
gated a lot in the past. The main problem with the sum on the 
right hand side is that it is dominated by the ultraviolet part of 
the spectrum and therefore is expected to have an ill defined 
continuum limit. But all spectral sums of the form 

s s (u) = E ^ E f( ZkX p) - E s f (tf; x ) 

k p x 

(9) 

k p=l 

define (nonlocal) order parameters for the center symmetry 
ifloll . Indeed, if we twist the gauge field with a center element 
z, we obtain: 

S f ('U) = zS f (U). (10) 

Our important observation is that, as the Polyakov loop, all 
spectral sums pick up a factor in the center of the group. Thus, 
not only the Polyakov loop, but also any other spectral sum 
of the above type might serve as a litmus paper for confine- 
ment. Of particular interest are sums which get their main 
contribution from the low lying eigenvalues. It has been con- 
vincingly demonstrated in iflQIl that the Gaussian sum with 
f(D) = cxp(— T>T> t) is very well suited for that purpose. For 
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a ST/ (3) lattice gauge theory, Figure[T]shows the Monte-Carlo 
averages of the partial sums 

3 n 

(?„((7)=^z^exp(-|^A p | 2 ) (11) 



k=l 



0=1 



and demonstrates their rapid convergence to a multiple of 
the (rotated) Polyakov loop expectation value. Actually it is 
sufficient to include less than 5% of the low lying eigenval- 
ues to obtain a decent approximation of the limiting value. 
Similar results hold for the spectral sums of the functions 
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FIG. 1 : Mean Gaussian sums 

QT for SU{3) on a 4 3 x 3-lattice near 
/3crit • The graphs are labeled with f3. 

f(T>) —V^ 1 and f(T>) =T>~ 2 corresponding to the propa- 
gators of T> and T> 2 . They are of particular interest since they 
relate to the celebrated Banks-Casher relation. 

The spectral problem (O is gauge-equivalent to the problem 
with twisted boundary conditions 



Vu z ip p = z X p z ip p with 

z ip p (x + N t , x) = -z* z ip p (x , x). 



(12) 



For calculating the gauge invariant spectral sums we may ei- 
ther twist the gauge field as in (O or twist the boundary condi- 
tions as in dT2b . Bilgici et al. extended the results in iflOTI and 
allowed for twists of the type z = exp(27rio;) in the bound- 
ary conditions 1 1311 . Admitting arbitrary values a £ [0, 1[, the 
twists are no longer center elements of SU (N) but only of 
U(N). Using nevertheless (fT2l as the primary definition of 
twist when z is extended to U(l) phases, the coefficients £„ 
of the expansion of the twisted quark propagator in powers of 
z can be easily obtained 111 311 . The coefficient £„ can be writ- 
ten as sums over all loops winding n-times around the torus 
in time direction. In particular, the coefficient Si is related to 
the spacetime integrated spectral density z g(X) for vanishing 
eigenvalue by 



Si = 



da 



z g(X \ 0), z = e 



2nia 



(13) 



In the numerical investigation below we will focus on trun- 
cations of local spectral sums (0. We will use IR sensitive 
spectral sums to find good approximations to (P(x)P(y)) or 
higher correlators. Guided by our previous results we consider 
the partial Gaussian sums 



g n (u ; x) = Y,4 ex P H Zfc M 

k p=l 



(14) 



where g p (x) has been introduced in (0 and (0. We will 
study by analytical (for cases which admit a complete analyt- 
ical evaluation of the spectral sums) and by numerical means 
whether and, in case, to which extent the correlator of two 
partial sums approaches the Polyakov loop correlator 



(Gn(U; x)G n (U; y)) — > const • (P(x)P(y)), 
when n tends to the total number n p of eigenvalues. 



(15) 



III. SPECTRAL SUMS AND CENTER SYMMETRY IN 
THE CONTINUUM 



So far the intriguing relations between spectral sums of 
twisted configurations and Polyakov loops have been estab- 
lished for lattice regulated gauge theories only. It is yet an 
open question which results remain meaningful in the contin- 
uum limit. Clearly, an object like Tt(V^') does not make 
sense in the continuum, and this was the main motivation in 
ifioll to introduce the generalized spectral sums (|9). Even if 
Tr f(D) exists in the continuum limit and even if, on the lat- 
tice, 



S f (U;x) = '£Sf(x ,x;U) 



(16) 



with Sf(x;U) from (0 is roughly proportional to the 
Polyakov loop, there are still the possibilities that this approx- 
imate proportionality is lost in the continuum limit or that the 
constant of proportionality diverges. For example, it was ob- 
served that the factor k in dU diverges in the continuum limit. 

In this section we study spectral sums for Euclidean gauge 
theories in the ab initio continuum formulation for a torus of 
extend /3 x L d_1 with L much bigger than the inverse tem- 
perature (3 = 1/ksT. The volume of the torus is V = (3 ■ V s 
while the spatial volume is given by V s = L d ~ x . On a torus 
the continuum Dirac operator 



(17) 



has discrete eigenvalues A p which are real for a vanishing 
quark mass. We shall consider Hermitian vector potentials 



A, = A^X a . 



(18) 



The gauge field is described by real- valued functions A® (x) 
(with a e {1, 2, . . . , dim(G)}) of the Euclidean space-time 
points x = (xq, x). The path integral measure contains fields 
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that obey periodic boundary conditions in the Euclidean time 
direction, 



A^Xq + f3, x) = A^xq, x). 



(19) 



Even for configurations with nonvanishing instanton number 
we may assume periodicity in the time direction lfl6[[l7ll . The 
Yang-Mills action is invariant under gauge transformations 



9 A^ = g(A„ + idjg- 1 , g(x) g G, 



(20) 



under which the field strength transforms as 9 F flL/ = 
gFfj, v g~ 1 ' In order to maintain the boundary condition $1% 
for the vector potential, we must demand that the gauge trans- 
formations are periodic up to a constant twist matrix z, 



g(x + P,x) = zg(x ,x). 



(21) 



When such a topologically nontrivial transformation is ap- 
plied to a strictly periodic vector potential A^ then 



(22) 



The gauge transformed potentials obey the boundary con- 
dition (fT9l only if z commutes with 9 A ll . This limits us 
to twist matrices in the center of the gauge group and ex- 
plains why twisted gauge transformations are called center 
transformations. In an irreducible and unitary representation 
of the group a center element is a multiple of the identity, 
z = phase factor ■ 1. The phase factor is such that z is a 
group element. We shall denote both the center element and 
its phase factor by z. 

In the absence of matter fields the twisted transformations 
form the global center symmetry. It can break spontaneously 
and the traced Polyakov loop 

P(x) = trV(x), V(x) = Vexp (^i drA (T,x)^, 

(23) 

which transforms nontrivially under center transformations 

9 P(x) = trPcxp (^ij^ ds 9 A (s,x)^j = zP(x) (24) 

serves as an order parameter for the center symmetry. As is 
well-known, the expectation value of P(x) is nonzero in the 
deconfining high-temperature phase and it is zero in the con- 
fining low-temperature phase. 

It is important to note that the center symmetry is explic- 
itly broken in the presence of matter fields in the fundamental 
representation. For example, quark fields transform as 



9 ip(x) = g(x)ip(x) 



(25) 



and are antiperiodic in Euclidean time. For a nontrivial twist 
the transformed field is not antiperiodic anymore, 



g ip(x +P,x) = -z' 1 9 ^(x , x). 



The eigenvalues of the Dirac operator for the transformed 
gauge field 



9 V A = V„ A = gV A g- x 



(27) 



are different to the eigenvalues of T>a if g is nonperiodic in 
time. Although 



V A ip p = Xpipp 



(28) 



implies 9 T>a 9 4'p = \ 9 '4 , pi f° r a nontrivial twist the 9 ifj p are 
no eigenmodes of 9 T>a because they are not antiperiodic in 
time. But for two gauge transformations g and g with the 
same twist z in (f2Tb the transformed operators do have identi- 
cal eigenvalues 



\ P { 9 V A ) = \ p ( 9 Va), 



(29) 



since the Dirac operators are gauge-related by the periodic 
gauge transformation gg^ 1 . 

Following the suggestion in 1 1 Qfl we consider the weighted 
sums 

S f (A) =J24 Tr/r ^) = J d d xS f {A-x) 

k 

S f {A;x)^Y,zt(x\trf^V A )\x) (3Q) 

OO 

k p=0 

but now for continuum Dirac operators. Here Tr denotes the 
trace over all degrees of freedom, whereas tr denotes the trace 
in spinor- and color space only. Similarly as on the lattice one 
collects the contribution to the spectral density of all eigen- 
function with the same energy, 



%{x) = Y,\ z i>vA*)Y 



(31) 



(26) 



The z q are gauge invariant scalars which transform nontriv- 
ially under center transformations. According to a theorem of 
H. Weyl lfl8tl the eigenvalues of T>a on a space of finite vol- 
ume have the asymptotic distribution \ p ~ p~ 1 / d such that 
the traces in (f30b exist for functions / which decay faster than 
l/\ d for large A. Actually, later we shall prove that the spec- 
tral sums defined as 

n 

S f (A;x)= Urn (32) 

p=0 k 

exist for a much bigger class of functions. 

The spectral function Sf transforms under center transfor- 
mations as follows, 

S f ( 9 A) = J24^f( Zk ^A) 

k 

= j24Trfr k v A ) (33) 

k 

= zY / 4'^f( Z ' k VA)=zS f (A), 
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where we have set zzk = zL The same argument applies to 
the density such that for all elements of the center we have 

S f ( a A) = z ■ S f (A) and S f ( 9 A; x)=z- S f (A; x). (34) 

All spectral sums Sf transform the same way as the Polyakov 
loop under center transformations and thus equally well serve 
as order parameters for the center symmetry. As on the lat- 
tice the eigenvalue problem for 9 T>a acting on antiperiodic 
functions is gauge-equivalent to 

V A ip p = Xptpp, ip(x + (3, x) = -z^ipfa, x). (35) 

For spectral problems the twisting of A^ has the same effect 
as twisting the boundary conditions with the inverse center 
element. Since the shifts \ p ( 9 T)a) — A p (2?a) only depend 
on the twist z we may as well choose a simple representative 
in every class of gauge transformations characterized by this 
twist. 

Gauge group SU(N): The cyclic center TL^ of this group 
is generated by 

z = cxp (2wiT/N) with T = diag (1, 1, . . . , 1, 1 — N). 

.(36) 

As simple gauge transformations with twist Zk ~- 
choose the powers h k of 

h(x ) = cxp (2irix T/PN) . 

The transformed gauge potential reads 



we 



(37) 



2irk 



Z «A^ = hK(xo)A„(x)h-"(xo) + — T5„, , 

k = l,...,N, (38) 
and the corresponding twisted Dirac operators are 



T> A = h k V A h~ 



(39) 



Gauge group U(N): The center of this group consists of 
the elements z = phase factor ■ 1 with arbitrary phase factors 
e 2ma . As simple representatives for the twisted gauge trans- 
formations with twist z we choose 

h a {x ) = exp(2TTiax /f3) ■ 1, < a < 1. (40) 

The gauge transformation shifts the potential by a constant 
proportional to the identity matrix, 



2ir 

z A^(x) = A^x) + —al 6^0, 



(41) 



similarly as in the construction of the Nahm-transform of self- 
dual U(N) -gauge fields lfl9tl . The twisted Dirac operators are 

"D A = h a V A h~ a , (42) 

and the sum over the center elements in the spectral function 
( f30b turns into an integral 



Sf{A) = J d d xSf(A- lX ), with 

S f (A;x)= f daz*(x\tvf( z V A )\x), 
Jo 



2-Kia 



(43) 

Although the center is not discrete the transformation rule for 
the spectral sums (f34-b applies. 



Comparing SU(N) and U(N): The equivalence of the 
spectral problems for z T>a on antiperiodic functions and T>a 
on functions with twisted boundary conditions (f35T > can be 
used to prove that the Dirac operators for certain su(N) and 
u(N) fields have identical spectra. To show this we consider a 
traceless potential A^ which can be viewed both as su(N) or 
as u(N) potential. We transform it with twisted gauge trans- 
formations g E SU(N) and g G U(N). The transformed 
potentials 9 A^ e su(N) and 9 A^ G u(N) are in general dif- 
ferent. However, if g and g are twisted with the same center 
element z of SU(N) then 



X P ( 9 V A ) = X P ( 9 V A ). 



(44) 



Note that for a = k/N the center element h a ((3) in d40b is 
actually in SU (N) and the result d44t applies. Thus for any 
su(A^)-potential A^ the Dirac operators with transformed po- 
tentials 



2nk 

+ J W TS^ esu(N) (45) 



A^(x) =A,(x) + ^lS^eu(N) 

have identical spectra. This observation is useful when one 
calculates spectral sums. 



A. Spectral sums and Polyakov Loop 

In the absence of matter a gauge invariant function is a func- 
tion of the gauge invariant Wilson loops based at some fixed 
base point, 



W Cx {A) = tr V exp i 



A^ 



(46) 



For a contractable loop these objects are invariant under cen- 
ter transformations and for loops winding A; -times around the 
periodic time directions they pick up the factor z k . We assume 
L 3> (3 in which case we may neglect Wilson loops winding 
around the spatial directions. 

Since Sf(A;x) is gauge invariant and transforms under 
center transformation the same way as a dressed _Pc„ with base 
point x = (0, £c), we conclude immediately that the functions 
F(C X , A) in the expansion 



Sf{A;x)=Y,Pe„-F{C x ,A) 



(47) 



are invariant under both twisted and periodic gauge transfor- 
mations. We conclude that for SU(N) these functions only 
depend on center symmetric Wilson loops based at (0, x). For 
example, they may still depend on Wilson loops winding A^- 
times around the periodic time direction. For the group U(N) 
with continuous center U(l) the function can only depend on 
contractable Wilson loops and center symmetric combinations 
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Actually, at least for the instanton solutions on the torus 
constructed by 't Hooft l20ll . a stronger result holds true, 
namely 



S f (A;x) ^ const -P(a;), 



(48) 



and this will be demonstrated in the following section. Since a 
similar relation approximately holds for certain spectral sums 
on the lattice one may conjecture that it holds for suitable cho- 
sen spectral sums in the continuum as well. 



B. On the convergence of spectral sums 

To investigate the convergence of the spectral sum we con- 
sider the Gaussian spectral sum 



Q'(t,A) = J d d xG'(t,A;x) 
g'(t,A;x) = Y,4(x\trcxp (-t z "V 2 A ) \x) - Q {x). 



(49) 

We subtracted the center-averaged density go (x) of the zero- 
modes for later use. More generally, if ij} p ^{x) are the or- 
thonormal eigenfunctions of T> A with eigenvalue /j, p = A 2 , 
then the center-averaged densities g p are 



k 



(50) 



with z g p (x) defined in QTb , wherein the z ip p ^(x) are eigen- 
functions of Z V A and not of z T>a as in the previous sections. 
In particular for gauge fields with a nonvanishing instanton 
number the zero-mode subtraction in d49l ) is always necessary 
and one deals with zero-mode subtracted heat kernels 



K'(t, A; x) = K(t, A; x) - g (x) = ^ e _AIpt 

p>o (51) 
with K(t, A; x) = (x\ tr exp (-tT> A ) \x) , 

where the sum extends over all p with p p > 0. On the torus the 
smallest nonvanishing eigenvalue \i\ is strictly positive and 
the zero-mode subtracted kernel falls off exponentially, 



K'(t,A;x) 



-tAI1 £>iO) for t 



(52) 



On the other hand, the heat kernel of the second order elliptic 



operator V A has the asymptotic small-t expansion (21 
1 N 

K(t, A;x) = j^{J2 + 0(t N+1 ) } , (53) 

where the Seeley-deWitt coefficients a n (x) are gauge invari- 
ant local functions built from the field strength and its covari- 
ant derivatives 12211 . Inserting this asymptotic expansion into 
the Mellin transform 



C(s,A;x) = wl dtt^K'^A^x) 
= E^p (s>d/2) 

p>0 



(54) 



one verifies that the analytic continuation C(s, A; x) of the 
sum in ( f54"b is a meromorphic function of s B23I1 . In even 
dimensions it has a finite number of simple poles as s = 
d/2, . . . , 2, 1 with residues ao(x), ai(x), . a^/2-i (x). Fur- 
ther its value at s = is 0^/2(2;) — Po(^) MM - Since the field 
strength and its covariant derivatives transform according to 
the adjoint transformation we conclude that the gauge invari- 
ant coefficients a n are invariant under center transformations, 
such that they cancel in the center average 

£' ( - 2s) ( x ) = YjsjJ dtt^G'it, A;x) = J2 4 C(«. Zk A, x) 

= £(E*i^-v4 (55) 



p>0 k 

It follows that E'( _2s ' (x) has actually no poles in the complex 
s-plane and that T,'^(x) = ~go(x). Without summing over 
the center elements the poles would not disappear and the sum 
overp would only exist for s > d/2. But if one first averages 
over the center and then sums over the quantum number p, 
then the last sum in (IB3T l exist/or all s. It is important that one 
first sums over the center elements and then over the p. For 
example, one finds 



d d a 



X">°\x) 



p>0 k 



d a x Zk g v 



» = 0, (56) 



since all densities z g p integrate to one. The same result follows 
from £''°'(.t) = — g~o( x ) since g>o integrates to zero. 

Let us now compare the spectral sums built from eigenval- 
ues and densities of T> 2 A with spectral sums built from those of 
T>a- An eigenvalue fi p of T> A is the square of an eigenvalue 
X p of T> A . In the massless case ±A p are both eigenvalues of 
the Dirac operator and using {T) A , 75} = it follows that the 
eigenvalue densities of T>a to positive and negative eigenval- 
ues are the same, g+ = g- = g where g is the density (l3TT l 
with eigenfunctions of T> A . Therefore we have g = 2g for 
the eigenvalue density of T>\. Thus we can conclude for the 
spectral sums of T>a 



t't- 2s \ X ) = £E^(r^r 2s + (- Zk \r 2s ) Zk e P 

p>0 k 

= ^EE^( 1 + (- 1 ) 2s ) \ zk v P r z %(x) 



p>0 k 



1 



(l + (-l) 2s )E'(- 2s) (x). 



(57) 



From this it is clear that in the massless case the spectral sum 
of T> A 2s exists in case the spectral sum of (V A )~ S exists and 
the spectral sums for V A S vanish for odd s. We conclude that 
also the spectral sums 



^) = E(E^t^^) 



p>0 k 



(58) 



where again the zero-mode contribution is omitted, exist for 
all s and gauge potentials A^. 
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IV. THE SCHWINGER MODEL AT FINITE 
TEMPERATURE 

For the Abelian instantons on the torus introduced by 
't Hooft ll20ll all eigenmodes of the operator V \ can be con- 
structed in the massless limit. The calculations for Abelian 
and non-Abelian gauge theories are very similar and so are 
the calculations in two and four dimensions. In this section 
we shall compute the spectral sums for all instanton configura- 
tions of the Schwinger model. We shall prove that the identity 
(|48l l holds true with a finite constant. 



A. Instantons and excited modes on the torus 

The [/(l)-gauge fields on the two-dimensional torus fall 
into topological sectors characterized by the instanton num- 
ber q. We choose a trivialization of the J7(l)-bundles such 
that in a given sector the fermionic field satisfies the 'bound- 
ary conditions' 



ip(x +(3, xi) = -ip(x ,xx), 
ijj(x , xi+L) = e^ipixo^x), 7 



2n_q (59) 



with q G Z. The gauge potentials are periodic up to a gauge 
transformation, 



A^xq^x + L) - A^xcXi) = d fl 'y(x). 



(60) 



The fields with minimal Euclidean action have constant field 
strength. We shall calculate spectral sums for the instanton 
solutions 



An 



<I> 

V 



Xi 



2nh $ 
— , Ai=0 with F 01 =B = —, 

(61) 

where V = (3L is the volume of spacetime and h an arbitrary 
constant. The flux $ of B is related to the instanton number 
by $ = 2nq. Without loss of generality we assume that the 
integer q is positive. The eigenvalues and q ground states of 
the massless Dirac o pera tor T>a = il^id^ — iA^) have been 
calculated earlier in 12511 . Here we shall construct all excited 
modes of T>a- 

The twisted gauge potential ( |4TT ) is equal to the untwisted 
potential with shifted h, 



a A = -Bx x 



2tt 
T 



"Aj = 0. 



(62) 



Hence it will do to study the spectral problem for A M in (f6Tb . 
The spectral sums will be compared with the Polyakov loop 
variable 



_ 2-jrih— i$ xi/L 



(63) 



For the instanton potential the straight and all dressed 
Polyakov loops with base (0, Xi) have the same value j63l . 
Then (|47| i implies that the spectral sums must have the form 



Below we shall calculate spectral sums for the squared 
Dirac operator T)\ with vanishing mass. For a positive in- 



stanton number its eigenvalues \i p = \ 2 are Il25ll 
(0 degeneracy: q 



Since they are independent of h they are also invariant un- 
der center transformations parameterized by a. It follows 
that Tr f( a T>) is independent of a such that the spectral sums 
Sf (A) in (l43l vanish for all functions /. Since the spatial 
average of P{x\) vanishes as well, this corroborates the con- 
jectured result d48l l. 

For the instantons ( f6Tb the spectral problem reads 



(65) 



with boundary conditions d59l . We choose a chiral represen- 
tation with 7°7 1 = (73 and diagonalize id^. On the cylinder 
[0, /?] x M the antiperiodic eigenmodes read 

XpA*) = e~ We 2 *"*"/" £ p ( Xl ), (£ G Z) (66) 

with time-independent mode functions £ p . These functions 
must solve the Schrodinger equation for the supersymmetric 
harmonic oscillator 

(-d 2 y + B 2 y 2 - Ba 3 ) £ p = 2 P B^ P , (67) 

where y is the shifted spatial coordinate 



y = x 1 + ^(£-h- 1/2) . 



(68) 



Eigenfunctions with <j^ p = £ p are called right-handed. Every 
right-handed eigenmode with energy 2pB gives rise to a left- 
handed eigenmode er 2 £p with energy (2p + 1)B. Hence we 
may as well focus on the right-handed sector. The zero energy 
states are 



B\ 



1/4 



-By 2 /2 



(69) 



where y(xi) has been defined in d68l ). The excited eigen- 
modes contain Hermite polynomials 



(70) 



X P A X ) = CpHpW B y) XoA x ), C P = 2p—[- 



Here we consider only right-handed modes and identify the 
nonvanishing component of a right-handed solution with the 
solution itself. The modes ( |70t do not satisfy the boundary 
conditions d59l since 



(71) 



X P A x o,xi +L) = e^^Xp^+gixo, 



but the true orthonormal eigenfunctions on the torus are just 
superpositions of these functions 



i>pA x ) = X P ,t+eq(x), £ = 1, 



,q. 



(72) 
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Note that the eigenvalues do not depend on I. Recalling the 
^-dependence of y the modes read 



xY, H p (y + sL ">) + aL) e 



a 2-nisqx / ' f3 



(73) 



We give a second representation which can be obtained by 
a Poisson resummation. In the appendix we show that the 
eigenmodes have the alternative representation 



P_ 2iri(l-l/2-qy/L)x /P 



^pA X ) 



(74) 

Forp = both sums define ^-functions and one recovers their 
modular transformation property. 



B. Spectral sums 

To calculate arbitrary spectral sum densities 
(x\ tr f(D\)\x) we determine the density of the eigen- 
value fi p in the right-handed sector, 



(75) 



The density in the left-handed sector is q~ = Q^_ v Hence- 
forth we shall skip the superscript +. To compute the sums 
over £ we use the representations ( f74b for the eigenmodes, 
since in this form they show a simple dependence on the quan- 
tum number £ (recall that y oc £). Using 



is proportional to P(xi). Due to the constraints imposed by 
the Kronecker symbol the double sum for g p , resulting from 
the series representation d74l i for the eigenmodes, reduces to 
one sum 



Q p (x) = c-a p (xo)- P(x\) with c = — 



1 IB 



L V 7T 



1/2 



, (79) 



where the series a p only depends on xq. Each term in the 
series contains a product of two Hermite polynomials H p 
and two Gaussian functions £o- Defining the variable x s = 
VB{xq + /3(l/2 + s/q)} the series takes the form 



a p (x ) = c 2 p e'^ 2 /' J2 H p {*> - ^/ 2 ) 

S 

x H p (x s + VBj3/2) e"< (80) 



To integrate over time we observe that a p is periodic in xo 
with period j3/q such that 



rP rP/l r 

/ dx j2--- = <i Yl — > 1 \ 

Jo Jo J- 



(81) 



Now we can apply the integral formula 1261 (7.377)] 



/oo nz 
dxH p {x + y)H p {x + z)e- x2 = \L p (-2zy), (82) 
-oo C p 

where L p denotes the Laguerre polynomial of order p, and 
this leads to 



dx a p (x ) = qj £ L p (BP 2 /2) er B ?l\ 



Inserting this result into ( F79b yields 



(83) 



e 2ni{s-r)y/L 

1=1 



= q5 s -r, nq e ^(s-r){ Xl /L-(h + l/2)/ q }^ n £ ^ (?6) 



the sum over I can be carried out. Twisting with a as in 
amounts to shifting h by a. Since the eigenvalues do not see 
the twist we may first calculate the integral of a g p over the 
center, 



g p {x) = / dae~^ la a g p (x) 



(77) 



and afterwards sum over p to calculate the spectral sums d43l . 
The a-dependence of a g p in d75l l comes only from the expo- 
nential factor in d76*| i and the corresponding integral over the 
center elements, 



dae -2-Kia e 2-Ki(s-r){x l /L-(h+a+l/2)/q} 



-P( Xl )S s - r ,- q (78) 



Qp{Xi) = / dx Q P (x) 



= -j-P(x 1 )L p (7rqr)e^/ 2 , r = |. 



(84) 



Taking the trace in spinor space amounts to adding the contri- 
butions of the right- and left-handed sectors. This finally leads 
to the following explicit result for the spectral sums in (l43l 

S f (A;x 1 ) = -^P(x 1 ) 

OO 

x fOh) { l pM + Lp-iinqr)} e^' 2 , (85) 

where we defined L_i = 0. This is the main result of this 
section. As expected on general grounds every function giv- 
ing rise to a convergent series (ISBI l defines a spectral function 
Sf(A, xi) which is proportional to the Polyakov loop. How 
fast the series converges to the asymptotic value depends on 
the particular choice of /. 
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Gaussian sum: Here we consider the Gaussian spectral 
sum 

G(A;xi) = J dx J dae~ 27Tla (x\ trcxp (- a V 2 A /pi 2 ) \x) 

(86) 

with some mass parameter pi. The integrand (x\ . . . \x) is just 
the heat kernel of a T> 2 on the diagonal in position space. The 
resulting series (l85l l with f(pi p ) = exp(— 2Bp/pi 2 ) relates 
to the generating function for the Laguerre polynomial I27I 
10.12(17)] 



T ^ex P (- I ^)=^L p (^ 

v ' p=0 

One obtains the simple result 

G(A; xi) 



\z\ < 1. 



(87) 



— coth — exp 
L " 2 



irqr 



/' 



coth — ) 



(88) 

Using B = 2irq/(3L one finds for L > qT/pi 2 that the re- 
lation between the Gaussian spectral sums and the Polyakov 
loop is the same in all instanton sectors, 

G(A;x 1 ) -ii^e-^'^PCari). (89) 

A natural energy scale at finite temperature would be the tem- 
perature itself, [i = T . With this choice the infinite-volume 
result simplifies further, 



Qoo{A\X\) 



T 



2 7 re 1 /4 



P{xx) for pi = T, 



(90) 

Propagator sum: Here we consider the propagator sum 

^~ 2 \A- Xl ) = J dx J dae~ 2ma (x\ tr' ( a V A y 2 \x} , 

(91) 

where tr' means the trace without singular contribution of the 
q zero-modes. Making use of the summation formulas 1261 
(8.976)], JH 



E 

P =i 



L p (x) 



P 



and Yl 

P =i 



L p -i{x) 



P 



-7 — logo; and y ] ^ p ^ y " J = e x T(0, x) 

(92) 

the spectral sum is given in terms of the Euler constant 7 and 
the incomplete Gamma-function, 

£(- 2 ) = A{ 7+ log( 7r(7 r)-e 7r9T r(0,7r g T)}e- 7r « T/2 F(.T 1 ). 

(93) 

In the large volume limit r = (3/L tends to zero and one 
obtains the simpler relation 

E( _ 2) £ lQg(7rgr)) . p{xi) + 0{(3/L y (94) 

In two dimensions Tr V^ 2 is logarithmically divergent in the 
ultraviolet for all background fields. For the instanton poten- 
tial with n p oc p this is evident. On the other hand, the spec- 
tral sum £( -2 ) is finite. Integrating over the center removes 
the divergence. 



O.i 



0.6 



y(-2) 



0.4 



0.2- 





g = 0.01 - 




q = 0.03 


T 


q = 0.1 





100 200 300 400 500 600 

n 

FIG. 2: Partial propagator sums e£ for P(xi) = 1 and their 
limiting values for different torus parameters (J3 = 1). 



Figure|2]shows the partial propagator sums 

^ { n 2) = J E - { L vM + Lp-xinqr)} e^' 2 (95) 

for P(xi) = 1 as a function of the included eigenmodes and 
their limiting values calculated from d93l ). Depending on the 
ratio r and the instanton number q they converge fast to their 
limiting values. 



V. FINITE TEMPERATURE SU(2) GAUGE THEORY 

In this section we calculate all eigenvalues and eigenfunc- 
tions of the Dirac operator for twisted and untwisted instan- 
ton configurations with constant field strength on the torus 
T 4 = [0,/3] x [0,L] 3 with volume V = (3 ■ V s . As a re- 
sult we obtain explicit expressions for the spectral densities. 
As expected, the Gaussian sums reproduce P(x) and get their 
main contribution from small eigenvalues. 



A. Instantons with constant field strength 

Following t'Hooft [20] we consider configurations with 
constant field strength, 



A 



-Ex 3 + ) 



(96) 

and assume that the constant chromo-electric and chromo- 
magnetic field components E and B are positive. The in- 
stanton number is proportional to EB and to the volume of 
space-time, 

1 f EB 

q = 32^ L e ""* tr F ^ = 2^ y ' (9?) 



J--1 
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and the Polyakov loop is periodic in the electric field, 

P(x) = 2 cos (2nh - Ef3x 3 ) . (98) 

For q > the gauge potential is periodic up to a nontrivial 
gauge transformation, A ll (x + L u e u ) = A^(x) + 8^"f u {x), 
with transition functions given by 

7i (.t) = BLx 2 cr 3 and 73 (x) = -ELx cr 3 - (99) 

The fermions are antiperiodic in time, periodic in X2 and fulfill 

i)(x + Le x ) = e^^ix), 



ip(x + Le 3 ) = e^^ipix). 



(100) 



Consistency demands that the fluxes in the 03 and 12 planes 
are both quantized, 

$03 = E(f3L) = 27rgo3, $12 = BL 2 = 2irq 12 (101) 

with (703, <7i2 G Z such that the instanton number q = 2<7o3<7i2 
is always even. In the chiral representation with 



ion 
I -mo I 



O, 



°o 



Oi J ' 7 \0 -<r J' 
(102) 

where era is the two-dimensional unit matrix, the squared 
Dirac operator takes the form 



^ = -D*-„ 38 (( B + E,<r ° ° ] 



with D 2 = D^D^ . The Pauli term acts with 0-3 on the color- 
S'[/(2)-indices, with (B + B)oo on right-handed spinors and 
with (B — E)(7o on left-handed spinors. D 2 is proportional to 
the four-dimensional identity in spinor space and commutes 
with do and 82- The (anti)periodic eigenfunctions decaying in 
the x 1 and x 3 directions have the form 

ip p/ (x) = e^ ixo / /} e 2m( - ioX °/' }+e2X2 / L ^ p (x 1 ,x 3 ). (104) 

On the functions £ p the operator —D 2 reduces to the sum 
of two commuting harmonic oscillator Schrodinger operators, 
one acting on x 3 and the other on x\. On an eigenfunction of 
03 in color space the operators read 



(-8 2 + E 2 y 2 )®l, 



H 12 = (~d 2 z +B 2 z 2 



i. 



(105) 



where we introduced the shifted coordinates (y, z) as follows, 

y = x 3 - — {h =F4 ± f } , z = xi T —^2 (106) 



903 



912 



for cr^p = ±£ P - Thus we recover two copies of the Schwinger 
model and conclude 

1> p j(x) = e -*./^(WW«/i)^( y )^ WXi 

(107) 



where \ is a constant right- or left-handed spinor and eigen- 
vector of er 3 in color space. The £'s are eigenfunctions of 
harmonic oscillator operators, 



#03^3 =2(p 3 + l)£Kp 3) 
Bi 2 £ P1 =2(pi + l)Se P i- 



(108) 



The eigenfunctions on the torus T 4 fulfilling the 'boundary 
conditions' dlOOl i are superpositions of the eigenmodes J 107b 
and read 



= '<P Pa ,£ ( X 0,X3)lpp 1 ,l 2 (x2,X 1 )x, 



(109) 



with £ 6 {1, . . . , 903} and ^2 G {1, . . . , (712}- The explicit 
form of the factors are given in d73l ) or ( f74b with obvious re- 
placements. For every pair of quantum numbers p = (p3,pi) 
there are qo 3 ■ ?12 eigenmodes of the squared Dirac operator. 

For cr 3 x — x in color space the variables y and z are given 
in j 106b with the upper signs and in spinor space the squared 
Dirac operator acts on the modes in ( 1109t as follows: 



V 2 A ^2(p 3 E + Pl B) I4 + 2 




Ea 



(HO) 



For CT3X = —x the variables y and z are given in (1 1061 ) with 
the lower signs and 



V\ -> 2 (p 3 B + jjiB + B) 1 4 + 2 



£0-0 




(HI) 



For every pair (pi , p 3 ) and generic E, B there are four eigen- 
values, each with degeneracy q = 2qo 3 • 912 ■ In particular, 
there exist q right-handed zero-modes in agreement with the 
index theorem. 



B. Spectral sums 

To calculate the densities g p (x) = J2e IVwWI 2 we use 
the representation (|74T i for the factors ipp3,e an( l ^piM m 
(1109b . Then the sums over £0 and £2 are calculated similarly 
as for the Schwinger model. The densities (x| tr f(T>A)\x) do 
not depend on X2 since T>a commutes with 82- Hence we may 
as well average over the ^-coordinate. 

This leads to a contribution 

Q P3 ( x 3)q P i(xi) = J dx Q P3 (xo,x 3 ) ■ i J dx 2 Qp 1 (x 1 ,x 2 ) 

(112) 

of the go3 • Q12 eigenmodes with fixed (p 3 ,pi) and fixed x- 
The explicit form of the factors is 

1 \ Q° 3 



00 ✓ 

(l + J2(-) n trV n (x)L P3 

n=l ^ 



Ef3 2 n 2 \ __ B ( / a„/ 2 ) 



0Pl \ X l) — j2 



\ + Y j tTQ n {x)L px 



BL 2 n 2 



-B(Ln/2) 2 



(113) 



I, (H4) 
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where we introduced the Wilson-loop winding once around 
the x 2 -direction, 

Q(x) =e l S MdX2 =&mg{e lBxiL ,e- iBxiL ). (115) 

Note that the factors do not depend on the choice of signs in 
( 1106b and this simplifies the analysis considerably. So far we 
did not sum or integrate over the twists and this is the reason 
why loops winding an arbitrary number of times around the 
time direction contribute to the sum d 1 1 3b . 

The eigenvalues of T) 2 A are given by the eigenvalues of the 
matrices in dl 101 111 II ). It is not difficult to see that eigen- 
functions with the following values of (p3,Pi) contribute to a 
given eigenvalue /i a .b = 2(aE + bB) of T>\ with a, b > 0: 

Ma.b (P3,Pi) = (a, b), (a - 1,6), (a, b— 1), (a — 1,6 — 1). 

(116) 

Generically there exist Aq eigenfunctions with the same eigen- 
value /j, a j,. But for a ~ or b = there exist only 2q eigen- 
functions and for a = and 6 = there exist q zero-modes. 
For a given eigenvalue /i a (, we have the following densities 

(a,6> 1) 

Qa,b{x) = 2{g a (x 3 ) + ga-i{x 3 )}{gb(xi) + g b -i(xi)} 

(117) 

where we defined p_i =0. 

Twisting in SU(2): Twisting the gauge potential inside 
the gauge group SU (2) amounts to adding 1/2 to ho or equiv- 
alently changing the sign of the Polyakov loop. The density 
Qb{x\) is unchanged but 



gQ3 

L 



Qa{x 3 ) - Z Qa(x 3 ) 

x £ trV«( X )L a (^)e- E 



ri=l,3,5, 



(n/3/2) 2 



(118) 



Only odd powers of the untraced Polyakov loop contribute to 
the sum g + z* z g over the center elements of SU(2). This 
confirms with our general analysis given earlier. 

Twisting in U(2): We may twist the su(2) gauge potential 
with center elements of U (2) or equivalently twist the bound- 
ary conditions by an arbitrary phase factor e 2rria . Averaging 
over the phases as in d77l i and below leads to 



J da a g a (x 3 ) e~ 2ma = P(x) L a ((3 2 E/2) e 



-E(0/2)\ 

(119) 

Integrating over the center of U(2) the sum over all windings 
n collapses to the contribution with just one winding which is 
proportional to P(x). Again this corroborates with our gen- 
eral analysis. 

Gaussian sum: To calculate the Gaussian spectral sum for 
the instanton configurations, 



Q(A; x)= [ dx Q 4 ( x \ trcx P {- Z ^A /M 2 
J k 



(120) 

one needs to calculate the sums ^ a . g a ^ exp(— fi a ^/n 2 ) 
with density dl 171 ). Summing over the SU (2)-center or inte- 
grating over the U (2) -center amounts to replacing g a in ( II 17b 



by the difference ( II 181 ) or the integral dl 19b . Again the result- 
ing series are calculated with the help of d87l ). After summing 
over the SU (2) -center elements one obtains 



q E B 

Q (A; x) = ttt- coth — coth ■ 



2V S y 2 



n=l,3,5,... 



£ (-rtr^Nexpj-^Lcoth-J 



E 



n' = l,2, 



1+ J2 trQ-'^expf-^f^coth^ 



(121) 



Not unexpected the first sum contains the center- and gauge 
invariant variables tr V , tr V 3 , tr V 5 . . . . If the spatial extend 
of the torus becomes large and we fix qo3 an d Q12, in which 
case the fields E and B tend to zero, then the Gaussian sum 
simplifies to 



M4«) = E 

n=l,3,. 



(-)™tr V n (x)-e-^ n/2)2 

(122) 



for L min{qT/n 2 , q/fj,}. On the other hand, for fixed 
q 03 , 3x2 ar, d fixed spatial extend we regain the the zero-mode 
contributions to (1121b for /.i~ 2 — > oo. For /i~ 2 — > we re- 
cover Qoo in ( 1122b - This implies an exponential decay with 
pt 2 as proposed in the general discussion on the convergence 
of the spectral sums. 

If we allow for [/(2)-twists with arbitrary phase factors 
then the resulting Gaussian sum is again given by the formula 
(1121b . but in the first sum over n only the term with n = 1 
contributes. In the thermodynamic limit L — > oo with fixed 
fluxes we find the simpler result 



g oo (A;x) = -£^^e-^ 2 ^P(x). (123) 



As expected, twisting with arbitrary center elements in U(2) 
removes the higher powers of the Polyakov loop. The formula 
is almost identical to the result d89l for the Schwinger model. 

Propagator sums The propagator sums are not absolutely 
convergent and the summation has to be carried out over fixed 
energy shells. Thus they are defined as 



Z'(~ 2s \x) 



lim 



,,s 

a,b 
0<Aiab<A 



Qa.b- 



(124) 



The existence of the right hand side for s > can explicitly 
be shown with the Mellin transformation of the zero mode 
subtracted heat-kernel (1121b in accordance with the general 
discussion following ( 
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FIG. 3: All eigenvalues of the staggered quark operator for a particular lattice configuration on a 6 4 lattice (left). Contribution to the interface 
energy as function of the mode index (right). 



TABLE I: Parameter set used for the simulations 



p 


7i 


72 


era' 2 


1.35 


2.0348 


-0.10121 


0.1244(7) 



where the phase factors are given by rj^x) = 



B. Energy of the Z2 interface 



VI. NUMERICAL INVESTIGATION 

A. Numerical setup 

Our lattice gauge simulations were carried out on a N t x 
Ng lattice using an improved action which is optimized for 
good rotational symmetry and good scaling ll29ll . We confined 
ourselves to the gauge group SU (2) and to a limited range of 
lattice spacings in this first exploratory study. The action is 
given by 



S 



(125) 



where P^ v {x) is the standard plaquette expressed in terms of 
the link fields U^x) € SU(2), i.e., 



1 



P^(x) = - tr U^x) U v (x + fi) Ul{x + v) Ut{x) 



(126) 



(2) 

and Pfiu (x) is half the trace of the 2 x 2 Wilson loop. We 



used the parameter set given in Table HI 12911 . Thereby, a is the 
lattice spacing and a is the string tension. For the study of 
the eigenmodes of the Dirac operator, we used the staggered 
Dirac operator: 



\ v u \y) = U^Sx+^y -Ul(x- n)S x 



The group SU (2) has only the two center elements 1 and 
z = — 1. The center transformation U — > Z U in (O does not 
change the (gluonic) action ( 11251 ). Defining the energy of the 
interface by the action difference, A — (S[ Z U] ~ S[U]), there 
is no penalty in the action for inserting such an interface. 

This situation changes when dynamical quarks are in- 
cluded: the quark determinant is not invariant under the map- 
ping U — > Z U. Considering the quark determinant as an inte- 
gral part of the total action, the energy of the interface is now 
given by 



A 



In 



det Z V, 



v 



det Vi 



(128) 



Representing each determinant det Vjj by the product of the 
eigenvalues A p of the corresponding Dirac operator, we can 
equally well write 



We note that in quenched approximation (quark effects on the 
link configurations are neglected), the surface energy A van- 
ishes since the configurations Z U and U contribute with equal 
probability to the Monte-Carlo average. Nevertheless it is in- 
structive to study a p for a single lattice configuration gener- 
ated in quenched approximation. Figure [3] left panel, shows 
all 2N t Ng eigenvalues for a particular lattice configuration, 
while the right panel shows the contribution a p to the inter- 
(127) face energy. We observe that the dominant contribution to 
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the interface energy arises from the low lying eigenmodes of 
the Dirac operator (Figure [3] right panel). Our findings also 
suggest that the mode sum in ( 1129b converges. This would 
imply that the interface energy is entirely determined by the 
IR regime of the quark sector. 



C. Polyakov loops 



Using the eigenvectors ip p (x) and eigenvalues X p of the 
quark operator, the Polyakov loop P(x) in (O can be recon- 
structed at position x by 



(130) 



P(x) = S n Jx ,x), Vx 



where n p — 2N t N^ is the total number of eigenmodes. It was 
already observed earlier that the mode sum in ( 1 1 301 > is domi- 
nated by the high end of the Dirac spectrum. Restricting the 
mode sum to a smaller upper limit n <C n p , we do not ex- 
pect that P{x) and S n {x) are correlated locally. In order to 
explore which value of n must be used to obtain a satisfactory 
agreement, we chose n = 0.9 n p and produced the scatter plot 
in FigurelU left panel. For n = n p , we observe a perfect cor- 
relation (which served as a benchmark test for our numerical 
approach). Already for n as high as 0.9 n p , this correlation 
has disappeared. A similar result holds for the expectation 
value of the Polyakov loop H §, [H . 

In order to quantify this correlation, we introduce 



P{x)S n (x , x) 



P*(x)) (S*(xo,x) 



(131) 



The average extends over the space-time index, and only con- 
tributions from a single lattice configuration are taken into ac- 
count. If P and <S„ are completely uncorrected, we find (in 
the confining phase) 



(P) 



0, 



0, uj n « 0. 



(132) 



By contrast, if both quantities are perfectly correlated, i.e., 
P oc S n , we obtain uj„ « 1. Figure[4] right panel, shows u> n 
as a function of n/n p , which is the fraction of the spectrum 
which was considered for the mode sum (1130b . Although the 
correlation is perfect for n/n p — > 1, a decent correlation is 
only achieved if almost all of the spectrum is taken into ac- 
count. 

Let us complete this subsection by replacing the mode sum 
S n in (1 1 301 > by the IR weighted sum 



Gn{x) :-- 



1 ™ 



z g p {x)e 



za2 / ; 



(133) 



In complete analogy to dl31t , we may define which quan- 
tifies the correlation between the Polyakov loop at each point 



TABLE II: simulation parameters 
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aa 2 


lattice 


T/T c 


configurations 


1.35 


0.1244(7) 


12 3 6 


0.7 


8658 


1.35 


0.1244(7) 


12 3 4 


1.0 


12000 


1.35 


0.1244(7) 


12 3 2 


2.1 


12000 



in space and G n (xo,x). Figure [3] shows our findings for 
/X = 1/3 and /i — > oo in comparison with u) n in d 1 3 1 b - The 
/x = 1/3 graph hardly shows a dependence on n/n p simply 
because the mode sum is already efficiently damped by the ex- 
ponential factor. In the intermediate range 0.4 < n/n p < 0.7, 
the mode sum utilizing // — > oo has a sizable correlation with 
the Polyakov loop. 



D. The static potential 

On the other hand, we do expect that the expectation value 
for the Polyakov loop correlator 



C(r) = (P{x)P(x + re 3 ) 



(134) 



is dominated by the IR regime of Yang-Mills theory at least 
for sufficiently large separations r. We therefore define the 
mode sum approximation to C(x, r) by 



C n( r ) = \Gn(x , X) Gn(x a , X + re 3 ) 



(135) 



Note that the expectation value on the right hand side of the 
latter equation does not depend on the particular choice for 
x due to translational invariance. From the above correlation 
functions, the static quark potential V(r) and its mode sum 
approximation V Q (r) is obtained from 

V(r)/T= -In C(r), V g (r)/T = - In C„ e (r), (136) 

where T is the temperature. In this first numerical study, we 
have chosen the exponential mode sum with fi 2 = 0.1 and 
truncated the mode sum by setting n = 50. Thus, the lowest 
50 eigenvalues contributed. 

We also investigated the deconfinement phase transition 
which takes place if the temperature exceeds the critical value 
T c ?s 0.69(2) y/a lL3(ill . Temperature is adjusted by varying 
the temporal extent N t of the N t x JV| lattice: 

(137) 



Our simulations parameters are summarized in Table iHl 

Figure |6l left panel, shows the potentials V Q (r) in com- 
parison with the full static potential V(r) in the confinement 
phase for Nt = 6. Most striking is that the potential V^(r) 
linearly rises at large distances r and therefore shows confin- 
ing behavior. At small distances r, the potential V s (r) is flat 
and does not show any sign of Coulomb law. This was not ex- 
pected since the 1 jr part arises from the exchange of gluons 
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m 





FIG. 4: Correlation between the Polyakov loop P(x) and the mode sum S n (x, 1) at each point a; of a single lattice configuration, 6 4 lattice, 
(left). The correlation measure lu„ dl31b as a function of n/n p , single lattice configuration, 6 4 lattice (right). 



o standard 
. u=l/3 
u=infinity 




-0.5 



FIG. 5: The correlation measure w n dl3U as a function of n/n p for 
different types of mode sums, single lattice configuration, 6 4 lattice 
(right). 



and belongs to the realm of the UV regime. The rather flat 
behavior implies a constant correlation at short distances and 
points toward a rather smooth field Q n {x). 

Also shown in Figure|6]are both potentials, i.e., V s (r) and 
V(r), in high temperature regime. At temperatures T m 2T C , 
both potential are essentially flat, and, in particular V s (r), has 
lost any signal of the linear rise. 

The situation is less clear for T w T c in Figure [6] right 



panel. While the Polyakov loop correlator can be, to a good 
extent, fitted by a 1/r Coulomb law, the mode sum still shows 
a significant linear rise. The interesting question is whether 
in the case of the mode sum the shift of the critical tempera- 
ture to higher values is nonvanishing in the continuum limit. 
To answer this question, much more time consuming simu- 
lations using higher values of (3 and therefore larger lattices 
are necessary. It would also be interesting to study the criti- 
cal temperature signalled by V s (r) for several types of mode 
sums. This is left to future work. 



E. Visualization 

The potential V Q (r) flattens for small quark antiquark dis- 
tances r. This indicates that the mode sums Q n (x) in (11331 ) are 
rather smooth functions of x. In order to get a first impression 
of the space texture of these mode sums, in Figure|7]we have 
visualized 



\Qn{x)l 



25, 



/i= 3 



within the spatial hypercube of a particular 20 3 x 8 lattice con- 
figuration. The lattice configuration was generated with the 
improved action JH] using (3 = 1.35. Using a = (440MeV) 2 
and V 'era 2 ~ 0.124, the length of the cube is roughly 3.2 fm. 
Thus, we observe a texture which is rather smooth at the 
length scale of 0.3 fm. This explains the behavior of the mode 
sum correlator ( 1135b at short distances. 
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FIG. 6: The static quark potential extracted from the truncated mode sum d 1 35 1 and from the Polyakov loop correlator i ll 34b in the confinement 
phase (N t = 6) and in deconfinement phase (Nt = 2) (left panel). The potentials for T ~ T c , (N t = 4) (right panel). 




FIG. 7: The mode sum m a 20 spatial hypercube. Side 

length of the cube: 3.2 fm. 



VII. CONCLUSIONS 

It is generally accepted that the low lying modes of the 
quark operator bear witness of the spontaneous breakdown of 
chiral symmetry in QCD. Here, we have investigated to which 
extent the low lying modes also contain information on quark 
confinement. 

For this purpose, Gattringer's mode sum approach @] to the 
Polyakov loop expectation value ^ P( x ) nas been general- 



ized to reconstruct the Polyakov loop locally. We have also 
pointed out that the mode sum approach is not restricted to 
lattice Dirac operators, but can be directly formulated in the 
continuum formulation of Yang-Mills theory. The existence 
of these sums has been studied in some detail. If one first sums 
over the center elements and afterwards over the eigenvalues, 
then the spectral sums exist/or all polynomial functions of the 
Dirac operator, especially for the IR dominated mode sums of 
interest, for example for 1/T>\. We have argued that the IR 
dominated mode sums equally well form an order parameter 
for confinement since these sums and the Polyakov loop share 
the same center transformation property. 

We have thoroughly investigated these mode sums by 
means of analytical calculations in the context of the 
Schwinger model and of SU(2) gauge theory in the back- 
ground of homogeneous field strength. As expected for these 
examples, the mode sums are proportional to the Polyakov 
loop for each point in space. 

Subsequently, we have employed SU (2) lattice gauge sim- 
ulations to study the relation between the low lying modes 
of the Dirac operator and the static quark antiquark potential. 
Below the critical deconfinement temperature, the correlator 
between two IR dominated mode sums is able to describe a 
linearly rising confining potential at large distances r. In the 
high temperature deconfined phase, this correlator reflects the 
flat behavior of the potential at large distances, and is in very 
good agreement with the correlator of two Polyakov lines. 
This clearly shows that the quark confinement mechanism is 
entirely encoded in the low lying spectrum of the Dirac op- 
erator. The search for confining degrees of freedom, such as 
vortices or monopoles, in the IR regime of the Dirac operator 
is left to future work. 



We finally point out that the preprint Bill appeared shortly 
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after we made this work available electronically. In this 
preprint, Bilgici and Gattringer also put forward the mode sum 
approach to the Polyakov line correlator and provided comple- 
mentary numerical insights for the gauge group SU(3). 
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and together with f p (rj) = £ P ( 7 1 + QV°) it leads to 

x £yaW f p ( m + w °) (A.9) 



For the q ground states the corresponding sums are Gaussian 
and give rise to theta functions 112511 . To perform the sums 
for the excited states we observe that under a Fourier trans- 
formation the step operators iA.3i are transformed into step 
operators, 



(at/)^) = iat/fa), 



(A. 10) 



where 



APPENDIX: POISSON RESUMMATION FOR HARMONIC 
OSCILLATOR EIGENFUNCTIONS 

In this appendix we consider and resum the eigenfunctions 
of the squared Dirac operator —T>\ for the Schwinger model 
on the torus. We measure Euclidean time in units of (3 and 
lengths in units of L. We assume that the flux $ = 2irq enter- 
ing the gauge potential 

A = -fei + 2-n-ho and A\ = 2nhi (A.l) 

is positive. The periodic eigenmodes on the cylinder [0, 1] x ffl 
have the form 



x P A x ) 



(A.2) 



with y 1 = x\ + i (£ — ho) , where I £ Z and the mode func- 
tions £ p are eigenfunctions of at a, 

a + a& = 2j>&, a=(2irq)iy 1 + (2irq)-id y i (A3) 

with [a, a'] = 2. The dependence on I enters via the l- 
dependence of y 1 . The normalized modes are 

MV) = c p a^Z (y) (A A) 

with 



1 



and eoW = (2g) 1/4 e 



1/4 p -Ttqx 2 



(A.5) 



p 2Pp\ 

The normalized eigenfunctions on the torus with boundary 
conditions (l59l l are superpositions of the modes on the cylin- 
der, 



e -~2iriqx xi 2iri(hoy —thi/q) 



(A.6) 



where we introduced the auxiliary function 

fp(v 1 ) = <P* n/>v %(v 1 ) with y° = Xo + ^-. (A.l) 
Here we can apply the Poisson resummation formula 



dy 1 e 2 ™™ 1 fiy 1 ), 



(A.8) 



a T = (2irq)i V - (2irq)-*d v (A. 11) 

with dual 'instanton number' q related to q by 

qq=l. (A. 12) 

The step operators a, at obey the same commutation relations 
as a, at. Since the ground state £o is transformed into the 
ground state with q we conclude that 



(A.13) 



with £o(v) — (2<?) 1 / 4 e 7T ^ V . In calculations it is advanta- 
geous to use Hermite polynomials H p generated by at and 



at 



a^ (y)=H p (y2^-yj^(y), 



In terms of these polynomials the equivalent series ( IA.6I ) take 
the form 



(A. 14) 



i/) Ptt (x) =c p e 2 * itx °e 2 * ihlXl 



x e 2mmqV °H P (y^Tq (y 1 + m)J £0(2/ + m) (A. 15) 

m 

and the resummed series JA.91 > reads 



1 C P e> -2mqxox lp 2m(hoy -ehi/q) 



- Ko (y ° + - ). 

3 7 9 

(A.16) 

In the main body of the paper we used antiperiodic eigen- 
modes of the Dirac operator. These are obtained by replacing 
h — > h + 1/2 and I — ► I — 1/2 in the results (IA.15I ) and 
(1X161) , 
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